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ABSTRACT 

Analytic modeling of several heat exchanger concepts to couple the SP-100 nuclear reactor primary lithium loop 
and the Space Stirling Power Convertor(SSPC)was performed. Four 25 kWe SSPC’s are used to produce the 
required 100 kW of electrical power. This design work focused on the inte rface between a single SSPC and the 
primary lithium loop. Manifolding to separate and collect the four channel flow was not modeled. This work 
modeled two separate types of heat exchanger interfaces (conductive coupling and radiative coupling) to explore 
their relative advantages and disadvantages. The minimum mass design of the conductively coupled concepts was 
18 kg or 0.73 kg/kWe for a single 25 kWe convertor. The minimum mass radiatively coupled concept was 41 kg or 
1.64 kg/kWe. The direct conduction heat exchanger provides a lighter weight system because of its ability to 
operate the Stirling convertor evaporator at higher heat fluxes than those attainable by the radiatively coupled 
systems. Additionally the conductively coupled concepts had relatively small volumes and provide potentially 
simpler assembly. Their disadvantages were the tight tolerances and material joining problems associated with this 
refractory to superalloy interface. The advantages of the radiatively coupled designs were the minimal material 
interface problems. 

INTRODUCTION 

During 1992 several concepts to integrate the SP-100 primary lithium loop to a Stirling convertor were explored. 
The SP-100 is a joint NASA, Department of Energy and Department of Defense space nuclear reactor system which 
passes lithium reactor coolant directly over the hot shoe of the multicouple thermoelectric power convertors. 
Stirling convertors have been proposed as a high efficiency power conversion alternative to the thermoelectrics. 
Originally a refractory alloy (1300 K) Stirling convertor was proposed for this purpose. Because of funding 
limitations the Stirling convertor development ended with a superalloy (Inconel 718) convertor. With a material 
substitution of Udimet 720 a 7 year convertor life could be achieved at a convertor hot end temperature of 1050 K. 
The Inconel convertor called the Space Stirling Power Convertor (SSPC) is being developed for the Civil Space 
Technology Initiative (CSTI) at NASA Lewis Research Center for this purpose. (Dudenhoefer 1992) The SSPC is 
designed to produce 25 kW of electrical power with an efficiency of 25% at a temperature ratio of 2.0 and a specific 
mass of 6.0 kg/kWe. 

This study focused on two separate types of heat exchangers to interface the SP-100 and the SSPC. The concepts 
are differentiated by the method in which the SP-100 loop and the Stirling hot end heat pipe condenser (called the 
Starfish Heater Head) axe connected. The first group of designs couple the hot primary lithium loop of the SP-100 to 
the Stirling convertor by radiating between the two systems. The second group couples the systems by a conductive 
interface similar to the design used for the multicouple thermoelectrics. Common to both are a lithium loop with 
refractory metal piping on the SP-100 side and a sodium filled, superalloy heat pipe on the Stirling convertor side. 
This heat pipe transfers heat to the working gas in the convertor. Several advantages and disadvantages for both 
ideas exist The primary advantage of the radiatively coupled designs Is that the separation of the heat source and 
the convertor prevented material compatibility problems. The primary problem with the radiatively coupled design 
is that the maximum heat flux achievable between the 1300 K lithium loop and the 1050 K Stirling convertor was 
about 8 W/cm^ while the evaporator of the sodium heat pipe has a heat flux limit of 40 W/cm^. This is due to the 
heat transfer limitations of Boltzman’s radiative heat transfer law for the temperatures being considered. 

A conductively coupled liquid lithium to Starfish heater head heat exchanger was explored to maximize heat 
transfer into the Stirling convertor. Although the conductively coupled designs have the potential of much greater 
heat fluxes, no single material is compatible with both liquid sodium on the convertor side and liquid lithium on the 


SP-100 side. Furthermore, heat transfer analysis shows that a heat flux in the neighborhood of 160 watts/cm 2 occurs 
when the hot primary lithium flows directly over an 1/32 inch Inconel heat pipe operating at 1050 K. This is much 
greater than the allowed 40 watt/cm 2 and must be compensated for in the design. 

This study developed two heat exchanger concepts for coupling the SP-100 to the SSPC. Calculations are made 
to predict heat flows, structural requirements, and thermal characteristics of the designs. Volume and mass estimates 
are made of the components for each heat exchanger. 

REQUIREMENTS 

Any of the heat exchanger designs proposed allow the Stirling convertor to be inserted directly into the flow path 
currently occupied by the thermoelectrics in the primary loop of the SP-100. Therefore, during the design work the 
heat exchanger configurations were made to match, as closely as possible, the current thermoelectric SP-100 
thermohydralic primary loop characteristics. These requirements set parameters such as temperature drop (75 K), 
reactor materials (Nb-l%Zr), and pressure loss (-2.0 psi). 

RADIATTVELY COUPLED DESIGN 

Analysis shows that radiatively coupled designs allow a heat flux of 8 watt/cm 2 into the evaporator of the 
Starfish while having the potential to operate at 40 watts/cm 2 . This dec rease in heat flux could have an adverse 
affect on both mass and volume of any of the radiatively coupled concepts considered. The primary advantage to 
the radiatively coupled design was that the separation of the heat source and the convertor prevent material 
compatibility problems. ' _~'Z ~ 

Several radiative concepts were developed in this design study. Of the designs considered the toroidal concept 
was found to be the lightest and easiest to manufacture. This paper will elaborate only on the toroidal concept This 
concept (See Figure 1 and 2), uses a single radiator panel made of Nb-l%Zr which is wrapped around the evaporator 
torus. The primary loop lithium flows through an inlet manifold that distributes the flow to Nb-l%Zr tubes 
embedded in the radiator panel which conducts the heat to the panel fins. The tubes are formed by the PLATECOIL 
fabrication process in the Nb-l%Zr radiator plate (Trantor Inc. 1990). In turn these fins radiate to the Stirling 
convertor’s heat pipe evaporators. The lithium is collected in an exit manifold and returned to the reactor. 

The evaporated sodium in the Starfish flows radially in to the convertor to the Starfish heater head and is 
condensed. The sodium is then pumped back to the evaporator through the capillary pumping action of arteries and 
then distributed to the evaporator through wire mesh wicking. Each sodium evaporator feeds half of the 25 kWe 
SSPC convertor. 

Thermal Analysis 

The SSPC helium gas temperature was required to have a hot side inner wall temperature of 1050 K. 
Calculations showed that a 16 K temperature drop can be expected from the inner wall of the Stirling convertor gas 
passages through the heat pipe to the outer wall of the evaporator. This required that the evaporator of the heat pipe 
operate at 1066 K. From the 1312 K average lithium temperature and the required 1066 K heat pipe evaporator 
temperature a total thermal resistance of 2.46 K/kW is calculated(Lubarsky 1955). 

An emissivity of 0.9 was used on both the super alloy Stirling convertor surface and the Nb-l%Zr side. These 
emissivities are the results of arc texturing research performed at NASA Lewis Research Center under the 
CSTI(Banks 1988). These values lead to a theoretical minimum required rad iator area of 1.28 m 2 to transfer lOO 
kWth between surfaces at an average of 1312 K on the lithium side and 1066 K on the super alloy side. All four 
concepts use a parallel tube configuration for the reactor side high temperature radiator panels. Each radiator panel 
lithium flow passage has an inner diameter of 5.08xl0" 3 meters. This results in Reynolds numbers of approximately 
10,000 and a total convective resistance of 1% of the allowable. Using a Nb-l%Zr radiator panel thickness of 
1.27xl0' 3 meters and a total of 20, 6.35x10' 3 meters diameter tubes, the total conductive resistance is estimated to 
be 10% of the allowable. With the addition of the radiati on resistance a radiator area of 1.5 m 2 is required . This is 
only 17% higher than the ideal radiator area of 1 .28 m 2 . 
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Multi-Foil Insulation 


Thermal insulation was needed to reduce the amount of heat lost into the convertor and to control the 
environment of the power system. Molybdenum Multi-Foil insulation was chosen because of its’ light weight and 
excellent insulation properties when used in vacuum. Estimates of heat loss through the insulation are made based 
on the temperature of the inner and outer layers of the Multi-Foil and the heat flux through the 
insulation(Worthaml988). A 5% heat loss set the required number of layers of insulation. The toroidal heat 
exchanger required 20 layers of Multi-Foil. 
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Stress Analysis 
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FIGURE 1. Cross Section of SSPC Radiatively Coupled Heat Exchanger. 


Calculations were made to evaluate required material thicknesses for the heat pipe in each of the concepts. The 
material used for the heat pipe is Inconel 617 that has an allowable stress of 62x1 0^ Pascal’s for 60,000 hrs of life. 
The maximum temperature this material is subjected to is 1066 K. To achieve the necessary area for heat transfer a 
0.308 meter torus was required. A burst disk was placed in the heat pipe wall to insure that if a high pressure helium 
leak occurs in the Stirling convertor/heat pipe interface the possibility of insertion of helium into the primary reactor 
loop is reduced. The burst disk was designed to operate at 1.48x10^ Pascal’s for the lifetime of the system without 
bursting. The internal pressure of the convertor is 1.55x10^ Pascal’s . The burst disk was sized to fail at 1.48x10^ 
Pascal’s . This lead to a 2.54x10'^ meter diameter disk. Mass results of the analysis are given in Table 1. Concept 
3 was chosen as the easiest to fabricate while also providing a lightweight interface between the SP-100 and the 
CIPC. 


Flow Channel Sizing 

Pressure drop calculations were made based on flow channel equivalent lengths around the reactor side of the 
liquid metal heat exchanger. Pressure drops included those in the inlet and outlet manifolds and the PLATECOIL 
heat exchangers. The lithium flow tubes were sized to give reasonably low flow velocities (< 3 m/s) to minimize 
pressure losses while always maintaining a turbulent flow (Re>4000) to enhance heat transfer. The wall thickness of 
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the tubes was set to accommodate the 0.151 Pascal’s pressure of the lithium loop. The selection of Nb-l%Zr as the 
lithium side material led to initial sizing of the tubing at 6.35x10”^ meters in diameter and 6.35x10"^ meters thick. 



FIGURE 2. Flow Path of Lithium in Radiatively Coupled Heat Exchanger 


Heat Pine 

A conservative 100 mesh stainless steel wire (0.01 inch‘d) was used for the wicking and arteries in all the 
concepts. This selection was based on preliminary results from wick erosion tests currently under test through a 
Lewis Research Center program. These results suggest that, without a coating, sodium erosion of the nickel in the 
stainless steel wicking is significant in thin wire screens. One hundred mesh wire was selected as the wicking size 
of choice because tests indicate this wire size will survive the seven year life requirement of the system. Because 
higher mesh sizes (smaller interwire spacing) provide increased capillary pumping ability it was desirable to use as 
fine a mesh wick as possible. Finer mesh sizes were considered and as more data becomes available this design 
choice may change. The 100 mesh wicking provides a total pumping height of about 0.254 meter’s in a lg 
environment . This was important because of ground testing considerations. 

Table 1. Weight Comparisons (Weight (kg) / Thickness (Meters)) 

Disks Radiator Fins Insulation Manifolds Total 
Single Torus 13.4 18.09 11.41 1.68 44 kg 

1.0E-3 1.27E-03 5.08E-04 8.9E-04 

In each of the designs wicking and artery attachment became a significant factor in the designs. Attachment 
methods for the wick and artery along with a relative ease of assembly influenced the choice of heat exchanger type. 
The single torus concept provided the easiest attachment of wicks and arteries . Because of the relatively large size 
of the torus in this concept a center separating plate was required between the two heat pipe sections to reduce the 
pumping height of the sodium when the system is ground tested. This effectively reduces the pumping height in 
half In all the concepts no additional mass penalty is included because of the wicking attachment 
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CONPUCnVELY COUPLED HEAT EXCHANGER 


A conductively coupled liquid lithium to Starfish heater head heat exchanger was modeled to increase heat 
transfer into the Stirling convertor and potentially reduce the heat exchanger mass. Although directly coupled 
Hpci gnc have the potential of much greater heat fluxes no single material is compatible with liquid sodium on the 
convertor side and liquid lithium on the SP-100 side. Furthermore, a heat transfer analysis showed that a heat flux 
of 160 watts/cm^ occurs when the hot primary lithium is passed directly over the evaporator of the Starfish heater 
head. This was much greater than the allowed 40 watt/cm^ and must be compensated for in the design. A unique 
solution to these problems is developed. Figure’s 3, 4, and 5 show the design in detail. Heat is transferred to each 
half of the SSPC separately by splitting the lithium flow in a manifold. Five cylindrical heat pipes evaporators (10 
total for the dual opposed convertor) on each side of the convertor are used to input heat to the heat pipe. The 
lithium is run in series through each of the five cylindrical heat pipes. Because the total temperature drop through 
the heat exchanger is 75 K, each heat pipe experiences a temperature drop of 15 K. 



Lithium 

Inlet 





Lithium 
Outlet 

FIGURE 4. Lithium How Channels for Heat Exchanger. 


Each of the proposed heat exchanger heat pipe evaporators consists of two basic components. The first is a 
Inconel 617 inner tube which acts as the evaporator section of the sodium Starfish heat pipe. To limit the heat flux 
into the heat pipes a thermal barrier coating (TBC) is applied to the outside of each of the Inconel tubes with a 
NiCrAly bond coat. Zirconia is used as the TBC because of its poor thermal conductivity and well-understood 
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material characteristics. This type of TBC has been used successfully for years on commercial jet engine turbines at 
much higher temperatures than will be seen in this application. Wicks and arteries are attached to the inside wall of 


Inconel 



FIGURE 5. Detailed Cross-Section of Heat Pipe Evaporator. 


the tube to distribute the sodium working fluid to the evaporative surface. The arteries run from the outside ends of 
the pipe to the roots of the fins in the Starfish to collect the condensed sodium and return it to the evaporator. 

The second component is the outer Nb-l%Zr shell which slip fits over the Inconel tube. This shell consists of 
two sub components. The first is the lithium flow passage that is comprised of two concentric cylinders that run the 
length of the inner Inconel heat pipe. At both ends of the concentric cylinders are donut shaped passages for the 
inlet and exit of the lithium. Thirty, equally spaced holes are made in the “donuts” at both ends to allow the lithium 
to enter and exit the concentric cylindrical lithium flow passages. This outer double shell is slip fit over the Inconel 
inner tube. Inconel’s coefficient of thermal expansion is much greater than that of Niobium and the two surfaces are 
in intimate contact when the system is brought up to operating temperature. Because lithium has a high affinity for 
oxygen and Nb-l%Zr is virtually transparent to the oxygen in the zirconia a method to slow the diffusion of oxygen 
through the Nb-l%Zr is needed. Molybdenum has been shown to slow the diffusion of oxygen and is used in the 
multicouple thermoelec tries on the SP-100 power system. A thin layer of molybdenum is Hot Isostatically Pressed 
(HIPed) to the Nb-l%Zr inner tube to slow this diffusion process. 

The need for the insulating layer of zirconia resulted from the maximum heat flux constraint of 40 W/cm 2 . To 
achieve the desired temperature drop of 75 K across the five series heat pipes it is necessary to have a minimum 
lithium inlet temperature of 1141 K (1066 K + 75 K). With only the Inconel and Nb-l%Zr layers present, the 
resulting heat flux into the first evaporator would be 130 W/cm 2 . Therefore, an insulating layer is needed to control 
the maximum heat flux into the evaporator. Since the lithium temperature entering each of the five series heat pipes 
is different (e.g., there is a 15 K temperature drop through each heat pipe) it is necessary to determine the zirconia 
thickness for each heat pipe that would result in the design goal of 35 W/cm 2 maximum heat flux. 

Thermal Analysis 

A computer model of a 5 degree section of the heat exchanger wall was constructed and used to perform the 
thermal analysis. PATRAN was used to generate the geometry and calculate the thermal resistance’s between 
thermal nodes. SINDA is used to solve for the steady state temperature solution. The lithium outlet temperature is 
calculated by the computer model, and the length of the heat pipe is varied until the lithium temperature drop 
through the heat pipe equaled 15 K (corresponding to 10 kW power input to each heat pipe). The computer model is 
run iteratively to determine the zirconia thicknesses, lithium temperature and heat pipe length that would result in a 
heat flux of 35 W/cm 2 at the entrance to each heat pipe and give the desired total heat input of 10 kW into each of 
the five series evaporators. A limit of 35 W/cm 2 was used because of the maximum tolerances achievable in the 
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spray coating of the zirconia. This limit can significantly alter the heat flux into the evaporator. A zirconia 
thickness of 2.03X10" 4 meters is used on the last of the five heat pipes as a minimum practical application thickness. 
The zirconia thickness for each of the four upstream heat pipes was then determined in sequence, using a lithium 
inlet temperature 15 K higher than that for the downstream neighbor. 

The main uncertainty in the analysis is the degree of thermal contact at the zirconia/molybdenum interface. This 
uncertainty was handled in the analysis by bracketing the thermal contact resistance between a conservative stainless 
steel/niobium interface value and no contact resistance. 

Stress Analysis 

This modeling and static analysis is done using PATRAN. The casing shell which is 0.3175 meter diameter and 
8.9X10" 4 meters thick is attached to the ribs. There are two end ribs and a center rib which are used to mount the 
casing shell to the Stirling convertor. These ribs are 2.54xl0‘ 2 meters deep radially and are an integral part of the 
convertor housing. The center rib is 3.2x10'^ meters thick while the end ribs are 4.7x10"^ meters thick. The center 
rib reduces the effective pumping height of the fluid which increases the heat pipe pumping margin. There are ten 
heat pipes located in two circumferential rows, equally spaced and staggered as shown in the figures. The heat pipes 
extensions are 8.89xl0' 2 meters in diameter, 0.13208 meters high from the casing shell and 8.9X10" 4 meters thick. 
The heat pipes are fitted with a heat source sleeve that has an initial diametral clearance of 3.5x10"^ meters between 
the heat pipe and the sleeve. This analysis reflects the state of stress with this initial clearance closed and an 
additional interference fit of 1.25x1 O' 4 meters between the Niobium sleeve and Inconel evaporator at operating 
temperatures. The heat pipe, casing shell, and ribs are made of Inconel 617. The casing shell, heat pipes, and ribs 
are modeled using QUAD elements. The loads are due to the 7.2720X10 4 Pascal’s internal sodium heat pipe 
operating pressure and the thermal coefficient of expansion mismatch between the Inconel and the niobium. The 
artnal temperature of the heat pipe and the sleeve is 1066 K and 1 150 K respectively,. The temperature stresses are 
due to any restrained thermal growth occurring after the clearance between the sleeve and the heat pipe is closed. 
The sleeve and heat pipe were modeled as if there is no clearance between them at a datum temperature that has no 
stress effect. To effect an interference fit of 1.262x1 0' 5 meters, the calculated required rise is temperature is 435 K. 
Because the temperature causes the entire convertor to expand uniformly it is fair to assume that the ribs are free to 
translate radially under the effect of temperature. For the pressure load, however, the housing, which is relatively 
stiff, may not expand proportionately. It may be fair, however to assume that the rib ends do not translate radially. 
Such situations offer difficulty in modeling the boundary conditions for the combined pressure and temperature load 
conditions. Since the thermal condition will have larger expansion than the pressure condition, the boundary 
conditions used were from thermal considerations. However a comparison of stresses due to this variation of 
boundary conditions were made. Table 2 shows a mass breakdown of the system. 

Table 2, Weight Comparisons (Weight (kg) / Thickness (Meters)) 



Torus/Heat Radiator Fins 

Insulation 

Manifolds Total 

Conductively 

Pipe 

7.4 

121 

8.18 

L5 24 

Coupled Heat 
Exchanger 

8.9E-04 

1.27E-03 

5.08E-04 

8.9E-04 


Comparison of the effect of the displacement boundary condition at the inner edges of the ribs showed a 
maximum stress of 3.8054xl0 7 Pascal’s . These stresses are at the shell to heat tube junction. It is possible that the 
shell or a portion of the heat tube can be thickened to reduce the stress at this junction. Such optimizations can be 
performed before finalizing the design. Analysis shows the stress values due to the temperature in the Niobium 
sleeve portion was 1.404xl0 8 Pascal’s. Combined loads of pressure and temperature are 1.405xl0 8 Pascal’s in the 
Niobium sleeve and 1.557xl0 7 Pascal’s in the Inconel 617 heat tube. The allowable stress for Inconel 617 is 
6.204xl0 8 Pascal’s at 1066 K temperature for a fatigue life of 60,000 hours. For Niobium the stress for 1 % creep in 
100,000 hours is 1.350xl0 7 Pascal’s at 1350 K. 
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Multi-Foil Insulation 

Prior to determining the zirconia thickness for each heat pipe it was necessary to determine the heat loss that 
could be expected through the Multi-Foil insulation. This was necessary since the heat loss would affect the required 
lithium temperature. This heat loss was calculated over the outer surface area of each heat pipe and compared to the 
desired heat load into the evaporator. 

Flow Channel Sizing 

Pressure drop calculations were performed to size the Niobium tubing diameter and wall thickness to match the 
primary loop SP-100 characteristics. The wall thickness was set by manufacturing considerations rather than the 
internal pressure of 1.516x10^ Pascal’s in the SP-100 lithium loop. To insert the heat exchanger directly into the 
SP-100 loop required a pressure drop of 2.689X10 4 Pascal’s . Calculations were made using equivalent hydraulic 
diameters and used to get initial estimates on pressure drops throughout the heat exchanger. A total pressure drop of 
2.716x10^ Pascal’s was found for the drops through the heat exchanger. The majority (over 1/2) of the pressure 
drop occurs at the 30 hole inlet to each of the heat pipe tubes. To reduce this pressure drop the lithium inlet tube 
position on each of the heat pipes can be changed. This will be reflected in future design work. 


Heat 

An analysis was performed to determine the wicking, artery, and flow channel sizes need to transfer the sodium 
from the condenser section of the Starfish to the evaporator surfaces and back to the condenser of the heat 
exchanger. One hundred mesh screen was chosen as the wick size of choice in this heat exchanger. Tests in a 
simulated 1/10 segment Starfish heat pipe for the CTPC has demonstrated performance to heat fluxes through the 
evaporator wicking in excess of 21 W/cm^. Experience of respected heat pipe vendors suggests that this could be 
extended to over 40 W/cm^* This value was used in all of the designs. Computations of the boiling limit with 100 
mesh screen verify this limit, albeit by theory that is not supported by much experimental data. At the same time, 
other computations show that the capillary pumping limit of two layers of 100 mesh screen is not reached. This is 
true even in ground testing where hydrostatic forces impose a penalty. Moreover, the screen should fill the wick full 
of fluid in the evaporator before evaporation causes liquid flow to commence. Additional computations showed that 
three layers of 140 mesh screen would be equally satisfactory, and perhaps better, with respect to the boiling and 
capillary pumping limits. However, until more is known about the corrosion life of high nickel alloys a wicking 
choice involving coarser screening is desirable, and dictates 100 mesh screen as the front runner in this design. 

In the meantime, it must be noted that the arteries of the present design use screen much finer than that in the 
wicking. This was to avoid the penetration of the artery wall by sodium vapor, which would cause artery 
malfunction. The resolution of this contradiction (wicking mesh versus artery mesh) depends upon a technology 
effort that is currently underway. This involves the complete coating of all surfaces of the high nickel alloy by a thin 
corrosion-resistant material currently being studied. Should this technique be successful, it is quite likely that 
updates of the present design would be able to use finer screen. Also, because of the very low artery pressure drops, 
it may be possible to make the larger artery diameters less. 

CONCLUSIONS 

Two types of heat exchangers which would couple the primary loop of the SP-100 and the heater head of the 
SSPC were explored. Masses for the radiatively coupled designs varied from 41 kg to 176 kg. Of the radiatively 
coupled concepts a single torus design was selected based on fabrication and assembly issues along with a mass 
analysis. The mass for this design is 42 kg or 1.64 kg/kWe. A single design was explored for the conductively 
coupled concept Significantly greater detail in the analysis was performed on this design and leads to a higher 
degree of accuracy in the mass estimated and design configuration. The conductively coupled concept weighed 18 
kg or 0.73 kg/kWe. 

The direct conduction concepts provide lighter weight systems because of the ability to run the Starfish 
evaporator at higher heat fluxes than those permitted by the radiatively coupled systems. The primary advantages of 
the radiatively coupled designs were the minimal material interface problems. The conductively coupled concepts, 
besides the reduction in mass, provide a relatively small volume along with an easier assembly process (especially in 
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wicking of the heat pipe). Their disadvantages are tight tolerances and material joining problems associated with the 
refractory to superalloy interface. 
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